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Abstract 

We study the impact of the fuh 0{as) QCD and 0{aew) electroweak (EW) radiative 
corrections to the e~^e~ — > tiZ^ process in the standard model (SM), and investigate the 
dependence of the lowest-order (LO), one-loop QCD and EW corrected cross sections on 
colliding energy ^/s. The LO, QCD and EW corrected spectrums of the invariant mass of 
tt-pair and the distributions of the transverse momenta of final top-quark and Z'^-boson are 
presented. The numerical results show that the one-loop QCD correction enhances the LO 
cross section, but the EW one-loop correction generally suppresses the LO cross section with 
our chosen parameters. In the case of mn = 120 GeV, the QCD relative corrections can 
reach 43.16% when ^/s = 500 GeV, while the EW relative corrections have the values of 
%, -4.07% and -5.42%, when = 500 GeV, 800 GeV, 1.2 TeV, respectively. 
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I. Introduction 



The standard model(SM)[T], [2] of elementary particle physics has provided a remarkably success- 
ful description of almost all available experimental data involving strong and electroweak (EW) 
interactions. While on the contrary, the Higgs-boson in the SM has not been discovered experi- 
mentally yet. That implies that the mechanism of the electroweak symmetry breaking (EWSB) 
giving masses to vector gauge bosons and fermions, remains a mystery. Moreover, the SM suf- 
fers from a few conceptional difficulties, such as hierarchy problem. This has promoted intense 
researches in building extension models in order to solve the hierarchy problem, and leaded to a 
rich phenomenology at present and future colliders. 

The top-quark was discovered by the CDF and DO collaborations at Fermilab Tevatron in 
1995 [3, 4J. It opens up a new research field of top physics, and confirms again the three-generation 
structure of the SM. But until now our knowledge about top-quark's properties has been still 
limited [5]. In particular, the couplings of the top-quark to the EW gauge bosons have not yet 
been directly measured in high precision. Since the top-quark is the heaviest particle discovered 
up to now[6l[7], it indicates that among all the observed interactions of elementary particles, the 
top-quark mass term breaks the EW gauge symmetry maximally, and the detailed physics of the 
top-quark may be significantly different from the predictions provided by the SM. In other words, 
the top-quark probably plays a special role in EWSB. The studies of top physics could illuminate 
the mechanism which breaks EW symmetry. Thus we may guess the new physics connected 
with EWSB could be firstly found in the precision measurements involving top-quark. A possible 
signature for new physics could be demonstrated in the deviation of the any of the tt7, ttZ and 
thW couplings from the predictions of the SM. Until now there have been many works which 
devote to the effects of new physics on the observables related to the top-quark couplings[8], |9]. 
All these studies indicate that the vector and axial form factors in the couplings of top-quarks 
and gauge boson could receive large corrections in certain models of dynamical EW breaking [lOj. 
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As we know that it is very hard to obtain the information about the ttV {V = 7) couphng 
from the precise measurement of the top-pair production e"'"e~ — ^ 7*/Z* — at a hnear colhder, 
because of the difficulty in distinguishing the contributions of tiZ^ and tt'-f couphngs. The probing 
of ttj couphng via 77 —>■ tt process was investigated in Ref . [1 1] . It is found that the 77 — * tt 
process is sensitive to tt'-f couphng. In references [121 [131 [SI, the QCD and EW corrections to 
top quark pair production via fusion of both polarized and unpolarized photons in the minimal 
supersymmetric standard model(MSSM), are presented. The corrections are found to be sizable. 
At hadron colliders, a measurement of the EW neutral couplings via qq ^ '~f* /Z* ti is hopeless 
due to the strong interaction process qq{gg) ^ g* ^ ti. Instead, they can be measured in 
QCD ttZ^/'-f production and radiative top quark decays in tt events {tt — > 'jW~^W~bb). Each of 
the processes is sensitive to the EW couplings between top-quark and the emitted Z*^-boson(or 
photon) [12] . The coupling of top-quark and Z^ gauge boson can be probed via the measurement 
of process e~^e~ — > tiZ^ at linear colliders. The process e'^e" ttZ^ at the tree-level has been 
already discussed in the literature [16j by Hagiwara in the context of the SM, while CP-violating 
effects in e~^e^ ttZ^ process were studied in the framework of the Two Higgs Doublets Model[T7] 
and with model independent effective Lagrangian|18j. Recently, the NLO QCD corrections to tiZ^ 
production at the LHC have been calculated in reference |19j . 

The International Linear Collider (ILC) is an efficient machine for precise experiments designed 
in e~^e~ colliding energy range of 200 GeV < y/s < 500 GeV, and it has an integrated luminosity 
of around 500 (fb)^^ in four years. It would be upgraded to ^/s ~ 1 TeV with an integrated 
luminosity of 1 {ab)~^ in three years [20j. This machine has sufficient energy to produce top-quarks, 
and be ideahy suited to precision studies of many top-quark properties due to the advantage of the 
measurement being carried out in a particularly clean environment. Thus the anomalous coupling 
between top-quarks and Z°-boson can be probed by measuring precisely the e~^e~ ttZ^ process 
at a linear collider. 

In this work we present the calculations involving the full one-loop QCD and EW corrections 
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Figure 1: The lowest order diagrams for the e+e —>■ ttZ^ process in the SM. 

to the process e"^e~ — > tiZ^ in the SM. The paper is arranged as follows: In Section II we give the 
analytical calculation description of the Born cross section. The calculations of full 0{as) QCD 
and 0{aew) electroweak radiative corrections to tiZ^ production at an e^e~ linear collider are 
provided in Section III and IV, respectively. In Section V we present some numerical results and 
discussions, and finally a short summary is given. 

II. LO calculation for e+e~ ttZ^ process 

We use the 't Hooft-Feynman gauge in the leading-order calculation, except when we verify the 
gauge invariance. The contribution to the cross section of process e+e~ tiZ^ in the SM at the 
lowest order is of the order 0{al^) with pure electroweak interactions. There are totally nine tree- 
level Feynman diagrams, which are drawn in Fig{TJ The Feynman diagrams in Fig{T] topologically 
can be divided into s-channel(Fig{T](l)-(5)) and t(u)-channel(Fig{T]^6)-(9)) diagrams. 
We define the notations for the process e^e~ — > ttZ^ as 

e+(pi) + e-{p2) ^ t{h) + i{k2) + Z\h), (2.1) 
where the four-momenta of incoming positron and electron are represented as pi and respec- 
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tively, and the four-momenta of outgoing top-quark, anti-top-quark and Z°-boson are denoted as 
ki, k2 and correspondingly. All these momenta obey the on-shell equations pf = P2 = f^l, 
kj = kl = mf and k^ = m|. 

The amplitudes of the corresponding tree-level Feynman diagrams for the process e~^e~ tiZ^ 
shown in FigUl are respectively expressed as 

= ^ " {s - ml) [ik!+ k^y - ^l^f (^0^(fc2)t;(pi)7. (1 - 4.^ - 75) 

X u{p2V{k,), (2.2) 

Qswcw s [{ki + ^3)2 -mi] ^ V 3 J 

X lMk2)v{pi)-i''u{p2)e^{k:,), (2.3) 

= " is - ml) [ik!+ k,Y - m]f^^'^^^ - - + + "^'^ 

X (^1 - - 75^ v{k2)v{p,)^, (1 - 4^2^ - 75) u{p2V{k,), (2.4) 

_2g3 / 8 

-^^^^ = 7 X ^ , .2 2T^(^i)7"(-^2-^3 + mt)7J 1- -s^-75 

X v{k2)v{pi)-i^u{p2)t''{k^), (2.5) 

X (1 - l^l^ - 75) t'(*2)il(Pi)7. (1 - isl- - 1-5) «(p2)e"(*3), (2.6) 

M^^^ = X ^,.2l( ^u{ki)Yv{k2)v{pi)Mh -h + me) 

GswCw [ki + k2)^ [{P2 - ks)^ - mil 



X 



^^{l-Asl.-^,)uip2)e^ik,), (2.7) 



X (l-4s^-75) (^2-1^3 + "^e)7M(l-4s^ -75) M(p2)e''(A;3), (2i 
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X 7,(l-44-75)Mb2)e^(A;3), (2.9) 

X 7^(1-4^2^-75) (-^ + 1^3 + me)7.(l-4s^- 75) M(p2)e''(/i;3). (2.10) 

Then the total amphtude at the lowest order can be obtained by summing up all the above 
amplitudes given in Eqs. fl2.2l) - fl2.10p . 
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A<o = ^7W«. (2.11) 



1=1 

The differential cross section for the process e+e~ tiZ^ at the tree-level with unpolarized 
incoming particles is then obtained as 

doo = ^Y. l-^orc^$3, (2.12) 

spins 

where the color number A^^ = 3, the summation is taken over the spins of initial and final particles. 
The factor ^ is due to taking average over the polarization states of the electron and positron. 
(i$3 is the three-particle phase space element defined as 



III. One-loop QCD corrections to the e^e ttZ^ process 

The QCD one-loop Feynman diagrams and corresponding amplitudes of the process e^e~ — tiZ^ 
are generated by using FeynArts3.2 [2T]. We present the box and triangle diagrams of one-loop 
QCD corrections to the e^e~ — tiZ^ process in FigJ2J We take the definitions of one-loop integral 
functions in Ref . [22] . and use mainly the FormCalc4.1 package [23] to calculate the amplitudes 
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of one-loop Feynman diagrams and get the results in a way well suited for further numerical or 
analytical evaluation. The calculation of the amplitude for one-loop diagram has been performed 
in the conventional 't Hooft-Feynman gauge by adopting dimension regularization scheme in which 
the dimensions of spinor and space-time manifolds are extended to D = 4 — 2e. The numerical 
calculations of the integral functions are implemented by using our in-house programs developed 
from the FF package[25- these programs, the numerical evaluations of one-point, two-point, 
three-point and four-point integrals are evaluated by using the expressions from Ref. |25j, and 
the scalar and tensor 5-point integrals are calculated exactly by using the approach presented in 
Ref. [26] ■ The phase space integration for the process e+e^ — > ttZ^ is implemented by adopting 
2to3.F subroutine in FormCalc4.1 package. In the calculation of the correction from the hard gluon 
emission process e^e^ — tiZ^g , we use CompHEP-4.4p3 programP?] to generate all the tree- 
level diagrams and their corresponding amplitudes automatically, and accomplish the four-body 
phase-space integration. In CompHEP-4.4p3 and 2to3.F programs the well-known adaptive Monte 
Carlo method (the VEGAS algorithm) is adopted for multi-particle phase space integration |28j. 
Then we can obtain the NLO QCD corrected total and differential cross sections. The amplitude 
of the process e~^e^ ttZ^ including virtual QCD corrections at 0{as) order can be expressed as 

Mqcd = M, + M^^cd. (3.1) 

where A^qc'd renormalized amplitude contributed by the QCD one-loop Feynman diagrams 

and their corresponding counterterms. The relevant renormalization constants for top field are 
defined as 

= (1 + \5Zi:^QCD)t\ to = (1 + \5Z^,QCD)t''- (3.2) 

Taking the on-mass-shell renormalized condition we get the 0{as) QCD contributed parts of the 
renormalization constants as 

dZ^^QCD = -ReEf^QCDimf) - nil^Re [^Iqcd{p^) + T^lQcoip') + 2SfQCD(p')] Ip-^?, (3.3) 
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where Re represents taking the real part of the loop integrals appearing in the self-energies only, 
and the renormalized top-quark irreducible QCD two-point function is defined as 

^t,QCDip') = t [i^PlS^qcd(p') + ^PB^f^Qcnip') + m,^QCDip')] (3.5) 

where a and (3 are the color indices of the top-quarks on the two sides of the self-energy diagram, 
and Pl,r = (1 =F 75)/2- The unrenormalized top-quark QCD self-energy parts at QCD one-loop 
order read 

^t,QCDip') = ^wcDip") = 1^ + [P, 0, m,] + 25i[p, 0, m,]) , (3.6) 

and 

^iQCDip") = 1^ (1 - 25o[p, 0, m,]) . (3.7) 
The corresponding contribution part to the cross section at 0{as) order can be written as 



' ' * spins 



(3. 



where A''^ = 3, Cp = 4/3, and the bar over summation recalls averaging over the spins of incoming 
particles. 

The virtual QCD correction contains ultraviolet (UV) and infrared (IR) divergences. The UV 
divergences from one-loop diagrams can be cancelled after the renormalization by adopting the 
dimensional regularization. We have verified the UV finiteness for the renormalized amplitude 
analytically and numerically. 

The IR divergences in the M.qcd originated from the contributions of virtual gluon ex- 
change in loops. These IR divergencies can be cancelled by the real gluon bremsstrahlung cor- 
rections in the soft gluon limit, due to the Kinoshita-Lee-Nauenberg (KLN) theorem [29j. In our 
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calculation these IR divergencies are regulated by a small gluon mass. The real gluon emission 
process is denoted as 

e+(pi) + e-(p2) -> t{h) + + Z°(A;3) + g{k), (3.9) 

where the real gluon radiates from top (anti-top) quark line. We adopt the phase-space-slicing 
method [31] to isolate the IR soft singularity. The cross section of the real gluon emission process 
(13. 9p is decomposed into soft and hard terms 

AaS^ = Aa2Jf + Aa^. (3.10) 

In practice, we take a very small value for the cutoff AEg/Ef, in the numerical calculations, and 
neglect the terms of order AEg/Eb in evaluating the following integration for soft contribution |22[ 



where Eg < AEg <^ y/s, AEg is the energy cutoff of the soft gluon, and the gluon energy can be 



obtained by Eg = y\k\'^ + m?g. 

The cancelation of the IR divergencies in the virtual contribution part and the soft gluon 
bremsstrahlung correction, can be verified during our calculation. Thus we get an IR finite cross 
section which is independent of the infinitesimal fictitious gluon mass m^. The hard gluon emission 
cross section Aa'^^J^ is calculated numerically by using Monte Carlo program in CompHEP-4.4p3 
package. 

The UV and IR finite total cross section of the process e+e^ — > tiZ^ including the 0{as) QCD 
contributions reads 

= ao + Aa«^^ = ao + Aa^^^ + Aaf.Jf + AaZ''^ = ^o(l + 5<^^^), (3.12) 

where Ao"*^*"^ and S'^'"^ are the QCD cross section correction and QCD relative correction at the 
order of C(as), respectively. 
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IV. One-loop EW corrections to the e+e ttZ^ process 

In this section, we present the calculation of the (9(ae^) order EW corrections to the process 
e^e^ ttZ^ . We use again the package FeynArts3.2 to generate the one-loop EW Feynman 
diagrams and the relevant amplitudes of the process e^e^ tiZ^ . The one-loop EW Feynman 
diagrams can be classified into self-energy, triangle, box, pentagon and counterterm diagrams. 
Some pentagon diagrams are depicted in FiglS] as a representative selection. In our EW correction 
calculation we adopt the t'Hooft-Feynman gauge and the same definitions of one-loop integral 
functions as in Ref . [22] . We take the dimensional regularization scheme to regularize the UV 
divergences in loop integrals, and assume that the CKM matrix is identity matrix and use the 
on-mass-shell(OMS) renormalization scheme ^2]- The relevant field and mass renormalization 
constants are defined as 

eo = (l + 5^e)e, TUf^o = nif + 6mf, m'^^z),o = ^w{z) + ^^w{z)^ Hq = {1 + ^6Zh)H, 
Ao = UzazZ + (1 + \^Zaa)A, Zo = (1 + Uzzz)Z + UzzaA, 



fi = (1 + bzf)f^, = (1 + Uzf)!^. (4.1) 



2 V 2 ' " ^ 2 ^ 2 

^5z^^)/^ = (1 + 

With the on-mass-shell conditions, we can obtain the related renormalized constants expressed as 



8mf = -^Re [E^(mJ) + (mj) + 22} 



,2 _ ^„^Wf2 \ r 2 



Sm'w = ReT^T i^tv)^ Smz = ReE§,^{ml), 6Zh = -Re — — 1^2=^2 



OZjaA — |p2=o, OZjzZ — — -Tie \p^=m\i 



SZzA = 2^i^, 5ZAz = -2Re^^^. (4.2) 



The charge renormalization constant and the counterterm of the parameter sw can be obtained 
by using following equations [22j: 



I Cw ^ Sw ^ V 
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We use Eqs.f l3.3l) - fl3.4p to calculate the renormalization constants of the fermion wave functions, 
but the top-quark QCD self-energies in these equations(S^Q(^^, T^fg^^ and Sfg^^) are replaced 
by the corresponding fermion EW self-energies(Sj^^y(^, ^f^EW ^f,Ew)^ respectively. And the 
explicit expressions of the relevant EW self-energies in the SM can be found in the Appendix B 
of Ref. |22j. The UV divergence appearing from the one-loop diagrams should be cancelled by the 
contributions of the counterterm diagrams. In our calculations, it was verified both analytically 
and numerically that the total cross sections including 0{aew) one-loop radiative corrections and 
the corresponding counterterm contributions, are UV finite. Then the one-loop level virtual EW 
corrections to the cross section can be expressed as 

= cro^r = / d$3 E (-MS-MvT) , (4.4) 

' ' ^ spins 

where pi is the c.m.s. three-momentum of the incoming positron, d$3 is the three-body phase 
space element, and the bar over summation denotes averaging over initial spins. J^^i^ is the 
renormalized total amplitude of all the one-loop EW level Feynman diagrams, including self- 
energy, vertex, box, pentagon and counterterm diagrams. 

Analogous to the calculation of the QCD correction as described in above section, we generate 
the EW one-loop Feynman diagrams and calculate their amplitudes for e~^e^ tiZ^ process by 
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using FeynArts3.2 and FormCalc4.1 packages. The numerical calculations of integral functions 
are performed by adopting our in-house library. The phase space of process e^e~ — > ttZ^ is 
integrated by using 2to3.F subroutine. In calculating the correction from hard photon emission 
process e+e" ttZ^j , the CompHEP-4.4p3 program is used to generate the tree-level diagrams 
and amplitudes, and carry out the integration in four-body phase-space. During the calculation 
of the EW corrections to the process e+e~ ttZ^ , the IR divergence originating from virtual 
photon correction should be canceled by the real photon bremsstrahlung correction in the soft 
photon hmit. We use also the phase-space-shcing method and divide the cross section of the 
real photon emission process, denoted as e+(pi)e~(p2) — ^ t{ki)t{k2)Z^{k3)^{k), into soft and hard 
parts. 

^(^real = ^^so/t + ^^feard = ^o((),o/t + Qhard)- (4.5 j 

By using the soft photon(_E^ < Aii^^) approximation, we get the contribution of the soft photon 
emission process expressed as 

Qtk\ Qtk2 Pi , P2 



(4.6) 



_ki ■ k k2 ■ k pi ■ k P2 ■ k_ 

in which AE^ is the energy cutoff of the soft photon and < AE-^ <^ y^, Qt = 2/3 is the 
electric charge of top quark, E^ — ■^\k\'^ + is the photon energy. Therefore, after integrating 
approximately over the soft photon phase space, we can obtain the analytical result of the soft 
photon emission corrections to e'^e~ ttZ^ . The cancellation of IR divergencies is verified and 
the results shows that the cross section for e+e" — > ttZ^g is independent on the infinitesimal 
photon mass m-y in our calculation. 

For the convenience in analyzing the origins of the EW radiation corrections, we split the full 
one-loop EW level correction into the QED correction part (Acr^) and the pure weak correction 
part (Acr'^). The QED correction part involves the contributions from the one-loop diagrams 
with virtual photon exchange in loop, the real photon emission process e+e~ — > ttZ^^ , and the 
pure photonic contribution of the renormalization constants of the related fermions. The rest of 
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the total EW corrections remains with the weak correction part, 
one-loop EW corrected total cross section as 



Then we can express the full 



(4.7) 



V. Numerical Results and Discussion 

In the following numerical calculation, we have used the 3-loop evolution of strong coupling con- 
stant as(/U^) in the MS scheme with parameters Aq{^ = 203.73 MeV, yielding a^^^{m\) = 
0.1176. We take the relevant parameters as^ [7] : aew(O)-^ = 137.036, me = 0.519991 MeV, 
= 105.6583692 MeV, = 1.77699 GeV, mw = 80.403 GeV, mz = 91.1876 GeV, rrit = 
172.5 GeV, rus = 150 GeV, = 1.2 GeV, nib = 4.7 GeV, rriu = ma = 66 MeV, sin^ Ow = 
1 — = 0.222549, and fi = rrit + \mz- There the light quark masses (m„ and m^) have 

the effective values which can reproduce the hadron contribution to the shift in the fine structure 
constant aetu(^|)[33]- With above parameter choice we get Q;s(/i^ = {rrit + \f^zY) = 0.1040. 
Since LEP II and the electroweak precision measurements have provided that the SM Higgs-boson 
exists in the mass range of 114.4 GeV < mu < 182 Gd^jSll ES], we include the relevant Higgs 
contributions with niH = 120 GeV as a representation. 

During our calculation of one-loop corrections, we have to fix the values of the fictitious masses 
of photon and gluon(as IR regulators and nig), and soft cutoff 5s = AEg(^^)/Eb except the input 
parameters mentioned above. In fact, the physical total cross section should be independent of 
these regulators and soft cutoff. We have verified the invariance of the cross section contributions 
at QCD and EW one-loop order, Acr^^f ^'^^ = Aa^ff''^^ + Aa^^-f within the calculation 
errors when the fictitious photon and gluon masses, and nig, vary from 10~^° GeV to 10~^ GeV 
in conditions of 5^ = 2 x lO^'^, = 120 GeV and y/s = 500 GeV. The relation between the 
0{as) QCD {0{aew) EW) correction and soft cutoff 6s demonstrates in Figs|l](a-b) (FigsIHt^a-b)), 
assuming nig = 10"^ GeV {m^ = 10"^ GeV), mn = 120 GeV and ^/s = 500 GeV. We can see 
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Figure 4: (a) The 0{as) QCD corrections to the cross section for e~^e~ — > tiZ^ as the functions of 
the soft cutoff 6s = AEg/Eb in conditions of rUg = 10"^ GeV, ttih = 120 GeV and = 500 GeV. 
(b) The amphfied curve marked with the calculation errors for Ao"*^*"^ of FigJl](a) versus 6s- 

in Fig|l](a) and Figl5]^a) that the curves for Acr^,/^ '^'^ + Acr^^^'^^ and Aa^^J^'^^ are strongly 
related to the soft cutoff 6s-, but both the total QCD and EW one- loop radiative corrections, 
^^QCD^EW^^ ^^QCD,STy ^ Aa^Jf '^^ + Aa£;?'^^), are independent of the cutoff 6s within 
the range of calculation errors as we expected. In FigSI^b) and FigJSI^b), we make the curves for 
/S^cr^CD Acr^^ greater in size and mark them with the calculation errors, respectively. In the 
further calculation, we fix the soft cutoff, fictitious photon and gluon masses as (5^ = 2 x 10~^, 
= 10-1 GeV and nig = IQ-^ GeV - 

In order to verify the reliability of our calculations for the tree-level cross section of process 
e"'"e~ tiZ^ , the QCD correction from hard gluon emission process e~''e~ — ^ tiZ^g and the 
EW correction from hard photon emission process e^e~ — > tiZ^'j , we evaluated their numerical 
results by adopting different gauges and software tools. In Tabled], we list these numerical results 
by taking itih = 120 GeV and y/s = 500 GeV, and using CompHEP-4.4p3 program[27] (in 
both Feynman and unitary gauges), FeynArts3.2/FormCalc4.1 [211 [23](in Feynman gauge, and 
the phase space integrations are implemented by using 2to3.F subroutine for e~^e~ ttZ^ process 
and in-house 2 — 4 phase space integration program for e^e^ tiZ^g and e^e^ — > tiZ^'y 
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Figure 5: (a) The 0{aeuj) EW corrections to the cross section for e^e~ tiZ^ as the functions of 
the soft cutoff Ss = AE^/Eb in conditions of = 10"^ GeV, mu = 120 GeV and ^/s = 500 GeV. 
(b) The amphfied curve marked with the calculation errors for Aa^^ of Figl5]^a) versus Ss- 

processes.) and Grace2.2.1[36j package(in Feynman gauge), separately. We can see the results 
are in mutual agreement within the phase space integration errors. We have calculated also the 
tree-level cross section for the e+e^ —>■ tiZ^ process by taking the same input parameters as in 
Ref.|16j. and got coincident numerical results. 

For the check of the calculations for one-loop diagrams, we used the LoopTools2.2 library and 
our in-house program(Both are in two combination cases respectively) to calculate the one-loop 
QCD and EW corrections to the process e+e" - ^ ttZ^ {Aa^^jf'^^ + Aa^f ^'^^) numerically, in 
conditions of mn = 120 GeV , ^/s = 500 GeV, 800 GeV, and other relevant parameters having the 
values mentioned above. The LoopTools2.2 is used in two ways: (I) with FF package(case I), (II) 
with routines adapted from A. Denner's bcanew.f(case II). The in-house programs are also applied 
in two ways: one is to use completely our created codes according to the expressions in Refs. [22l[26] 
for the numerical calculations of N-point(A^ = 1, 2, 3, 4, 5) integrals (case I), another way is to use 
our in-house program according to the formulas in Ref. [26j for the numerical calculations of 5-point 
integrals based on FF package(case II). In all cases we applied 2to3.F subroutine in FormCalc4.1 
package to perform the phase space integrations. The numerical results are listed in Table [2j Since 
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CompHEP 


CompHEP 


Fey n Arts 


Grace 




Feynman Gauge 


Unitary Gauge 


Feynman Gauge 


Feynman Gauge 


ao (fb) 


1.0649(1) 


1.0648(1) 


1.0648(2) 


1.0647(2) 


A^.'i';." (fb) 


7.028(2) X 10-2 


7.028(2) X 10-2 


7.028(2) X 10-2 


7.028(2) X 10-2 


Aa^r, (fb) 


0.4631(2) 


0.4632(2) 


0.4631(2) 


0.4632(2) 



Table 1: The comparison of the results for the tree- level cross section for the 
process e~^e~ —>■ tiZ^ (o"o) and the corrections from hard gluon/photon emission 
processes e+e" — > tiZ^g {^cr^aiS) and e~^e~' ttZ^'-f {^(7^^^), in conditions of 
vriH = 120 GeV , ^/s = 500 GeV and 6s = 2 x 10"^. The numerical results are 
obtained by using CompHEP-4.4p3(in both Feynman and unitary gauges), Fey- 
nArts3.2/FormCalc4.1(in Feynman gauge, the phase space integrations are imple- 
mented by using 2to3.F subroutine for e~^e~ — » tiZ^ process and in-house 2 — 4 
phase space integration program for e~^e~ —>■ tiZ^g and e~^e~ — > tiZ^'j processes) 
and Grace2.2.1(in Feynman gauge) packages, separately. 

among all the one-loop QCD Feynman diagrams for e~^e~ —>■ tiZ^ process there is no pentagon 
diagram, we get the same results when we use both in-house program of case (II) and LoopTools2.2 
of case (I). Table [2] shows there exists the coincidence between the corresponding results within 
the calculation errors. 

In Figl6](a) we present the LO and 0{as) QCD corrected cross sections(o"o, 
Acr*^*^^) as the functions of colliding energy ^/s with rriH = 120 GeV. The corresponding relative 
QCD radiative correction(5*^'"'^ = ^^^^) is presented in Figl6]^b). We can see from Figsj6]J^a-b) 
that both the LO and QCD corrected cross sections are sensitive to the colliding energy when ^/s 
is less than 700 GeV, while decrease slowly when ^/s > 900 GeV. Figj6t^b) shows that the QCD 
relative radiative correction has a large value in the vicinity where the colliding energy is close to 
the tiZ^ threshold due to Coulomb singularity effect. 

In FigsJ7](a-b) we depict the curves of the LO cross section(cro), 0{aew) QED, weak, and total 
EW corrected cross sections((T^''^''^^ = ctq + Aa^'^'^^) as the functions of colliding energy ^/s 
with niH = 120 GeV, and the corresponding relative radiative corrections(5'^'^'^^ = Aof^^f!!!) 
versus colliding energy are drawn in Figd^b). FigsI7](a-b) show that the QED and total EW 
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LoopTools2.2(I) 


LoopTools2.2(II) 


in-house(I) 


in-house(II) 


vir+soft 


500 


0.3893(3) 


0.3894(3) 


0.3892(3) 


0.3893(3) 


'-^^vir+soft 


800 


-1.396(1) 


-1.396(1) 


-1.396(1) 


-1.396(1) 


— ^" uir+sofL 


500 


-0.5550(5) 


-0.5551(5) 


-0.5550(5) 


-0.5551(5) 


AfjEW 

vir+soft 


800 


-2.996(3) 


-2.996(3) 


-2.995(3) 


-2.996(3) 



Table 2: The comparison of the results for the QCD and EW one-loop corrections 
to the process e+e" ^ ttZ^ (A^?o/f + ^(^T^ and AagJ^ + Aa^^^ ) numerically, 
in conditions of ttlh = 120 GeV, ^/s = 500 GeV, 800 GeV, and other relevant 
parameters having the values mentioned above. The numerical results are obtained 
by using LoopTools2.2 library in two combinations ((I) with FF package, (II) with 
routines adapted from A. Denner's bcanew.f), our developed in-house codes for all 
the numerical calculations of N-point(A^ = 1,2,3,4,5) integrals (case I), and our 
developed program for the numerical calculations of five-point integrals based on 
FF package(case II). 




Figure 6: (a) The LO and QCD corrected cross sections for the process e'^e" — > ttZ^ as the 
functions of colliding energy ^/s with itih — 120 GeV. (b) The corresponding relative QCD 
radiative correction (5*^*-^^) versus a/s. 
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(GeV) 










5^''^(%) 


500 


1.0648(2) 


1.5244(4) 


43.16(3) 


0.9729(6) 


-8.63(5) 


600 


2.9432(4) 


3.5741(7) 


21.46(2) 


2.805(2) 


-4.70(7) 


800 


4.243(1) 


4.698(2) 


10.72(3) 


4.070(4) 


-4.07(7) 


1000 


4.187(1) 


4.478(2) 


6.93(4) 


3.994(4) 


-4.63(8) 


1200 


3.811(1) 


4.001(2) 


4.99(5) 


3.604(4) 


-5.42(8) 



Table 3: The numerical results of the tree-level, one-loop QCD, EW corrected cross 
sections and the relative QCD, EW radiative correct ions ((5*5C^.-^w^ = Aa'^'^^'^^ /ao) 
for the process e~^e~ ttZ^ , by taking niH = 120 GeV and ^/s = 500 GeV, 
600 GeV, 800 GeV, 1000 GeV, 1200 GeV, separately. 

corrections always suppress the LO cross section, while the weak correction part enhances the LO 
cross section, except it slightly suppresses the LO cross section in the region of ^/s > 1020 GeV. 
The QED relative correction is about —38.0% at tiZ^ threshold and rises to —2.8% at 1.2 TeV, 
but the weak relative correction is about 17.6% at threshold and decreases to —2.6% at 1.2 TeV. 
Thus the QED and weak contributions partially compensate each other when they combine to 
form the total EW correction, and yield the EW relative corrections of about —20.4% at threshold 
and —5.4% at 1.2 TeV. Again, we see the Coulomb singularity effects on the curves of corrections 
6^ and in FigjTl^b) in the vicinity of threshold. There the curves for the QED, weak, and total 
EW relative corrections in Fig|7|(b) demonstrate that the large EW corrections near the threshold 
region are mainly contributed by QED corrections which originate from the diagrams involving 
an instantaneous virtual photon in loop with a small spatial momentum. To show the numerical 
results presented in Figs|6t^a-b) and FigsjTl^a-b) more precisely, we list some typical numerical 
results for Born cross section(cro), QCD, EW corrected cross sections(cr*5'"'^, a^^) and QCD, EW 
relative corrections(5'^'"^'^'^ = Aa^'"^'^^ /ctq) for the process e~^e~ — > tiZ^ in Table [31 All these 
results show that when y/s goes up from 500 GeV to 1.2 TeV, the relative QCD(EW) radiative 
correction, 5^^^{5^^), varies from 43.16%(-8.63%) to 4.99%(-5.42%). 

In Figj8](a) we present the spectrums of the invariant mass of top-quark pair(Mtf) at LO, 
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500 600 700 800 900 1000 1100 1200 
Vs [GeV] 



500 600 700 800 900 1000 1100 1200 
Vs [GeV] 



Figure 7: (a) The LO cross section, QED, weak, and total EW corrected cross sections(cr^''^''^^) 
for the process e~^e~ — > ttZ^ as the functions of colhding energy ^/s with mfj = 120 GeV. (b) 
The corresponding relative QED, weak, and total EW radiative corrections ((J-^'^'^^) versus ^/s. 

QCD(EW) one-loop order for the process e'^e~ — tiZ^ by taking rriH = 120 GeV and ^/s = 
500 GeV. The distributions of the transverse momenta of top-quark((i(To/(ip^, da^'"^'^^ /dp^rp) and 
Z^-hoson{dao/dp§,, da^^^'^^ /dp^) are depicted in FiglHl^b) and (c), separately. Fig|8]^a) shows 
that the QCD corrections always significantly enhance the LO differential cross section dao/dMu 
in the whole plotted region, except the enhancement becomes much smaller when Mu > 400 GeV, 
while the EW correction either slightly enhances the LO differential cross section of dao/dM^ 
in the region of Mu < 365 GeV, or distinctly suppresses the differential cross section when 
Mft > 375 GeV as shown in FigJHJ^a). All these three figures demonstrate that the QCD corrections 
induce the enhancement to the LO differential cross sections of the invariant mass of top-quark 
pair, the transverse momenta of top-quark and Z'^-boson. We can see from Figsl8](b) and (c) that 
the EW radiative corrections can only slightly reduce the LO differential cross sections of da^/dptp 
and da^/dp^ in all possible and regions. 
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Figure 8: The spectrums for e+e" —>■ tiZ^ process at LO, QCD and EW one-loop order in 
conditions of mn — 120 GeV and ^/s — 500 GeV . (a) The spectrums of invariant mass of 
top-quark pair M^f. (b) The spectrums of the transverse momentum of top-quark p\,. (c) The 
spectrums of the transverse momentum of Z°-boson pf. 

VI. Summary 

Probing precisely the top physics and searching for the signature of new physics are important 
tasks in present and future high energy physics experiments. The future e^e~ linear collider could 
provide much more efficient laboratory to put these measurements into practice with a cleanest 
environment. In this paper we have studied the full one-loop QCD and EW corrections to the 
e+e~ — > ttZ^ process in the SM. We investigate the dependence of the effects coming from the 
QCD and EW contributions to the cross section of process e+e~ — > on colliding energy ^/s. 
We find that the QCD correction enhances the tree-level cross section, while the one-loop EW 
correction suppresses the LO cross section. Our numerical results show that when mu — 120 GeV 
and colliding energy has the values of 500 GeV and 1.2 TeV, the corresponding relative QCD(EW) 
corrections are 43.16%(-8.63%) and 4.99%(-5.42%), respectively. We present also the LO, QCD 
and EW corrected differential cross sections of transverse momenta of final top-quark and Z^- 
boson, and the distributions of top-quark pair invariant mass. We conclude that both the QCD 
and EW radiative corrections have relevant impact on the e+e~ — > tiZ^ process, and should be 
included in any reliable analysis. 
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